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We point out that by comparing the total mass (in gravitational units) of the earth- 
moon system, as determined by lunar laser ranging, with the sum of the lunar mass as 
independently determined by its gravitational action on satellites or asteroids, and the earth 
mass, as determined by the LAGEOS geodetic survey satellite, one can get a direct measure 
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of the mass of earth-bound dark matter lying between the radius of the moon's orbit and the 



must be less than 4 x 10 9 of the earth's mass. 



O j. geodetic satellite orbit. Current data show that the mass of such earth-bound dark matter 

Current interest in dark matter has been heightened by the recent report by the DAMA/LIBRA 
Q^i' collaboration [l[ of evidence for galactic halo dark matter, based on their observation of an annual 

2 ■ modulation signal. Astrophysical arguments suggest that the galactic halo dark matter mass 

density is around 0.3(GeV/c 2 )cm~ 3 , but it is still an open question whether in addition to dark 
matter bound to the galaxy, there may be larger dark matter concentrations bound to the sun, and 
bound to the earth. T : 
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le possibility of sun-bound dark matter was discussed in an article of Frere, 
who pointed out that local dark matter concentrations in the galaxy may 



OO 

have played a role in the formation of the solar system. Their paper, and the papers of Sereno and 
g ^ I of forio Q, and of K„ and P itj eva fl, ose _ t s b a S ed oo p,^ o rbitS 
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to place a limit on a local excess of sun-bound dark matter of order 3 x 10° times the galactic 
halo dark matter mass density. It is also possible that there may be further local concentrations 
of earth-bound dark matter, which if large enough could be relevant for understanding the 
recently reported [7|] spacecraft flyby anomaly. Thus, it would be useful to have a direct method 
for measuring, or at least placing limits on, the mass of earth-bound dark matter. 

The aim of this note is to show that it is possible to set a direct limit on the total earth-bound 
dark matter mass lying between the radius of the moon's orbit and the radius of low lying satellite 
orbits, such as that of the accurately monitored [s| LAGEOS satellite. For a satellite of negligible 
mass in a circular orbit around an astronomical body of mass M, measurement of the orbit radius 
R and the orbital period T gives the product GM (with G the Newton gravitation constant) by 
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use of the formula 

4vr 2 i? 3 

GM = - f ^- . (1) 

Thus, from a measurement of the radius and period of the LAGEOS orbit, one gets GM©, where 
we have defined the earth mass M© here to include the mass of all earth-bound dark matter 
lying within the radius of the LAGEOS orbit. Similarly, by a measurement of the orbit and 
period of lunar orbiters close to the moon [{J, one gets GM m , where here we have defined the 
lunar mass M m to include the mass of all moon-bound dark matter lying within the lunar orbiter 
radius; we shall assume this moon-bound dark matter mass to be negligible for purposes of this 
analysis. An alternative way of independently determining the moon's mass is to study the orbit 
of a near-passing asteroid, such as Eros [10], which is influenced by the gravitational field of the 
moon as well as the earth. From such an analysis one can extract an accurate figure for the 
ratio R & / m = (GM© + GAM e )/(GM m + GAM m ), with AM© and AM m denoting respectively 
possible small contributions from earth-bound and moon-bound dark matter. Again assuming that 
moon-bound dark matter can be neglected, and expanding in the small quantity GAM©, this ratio 
becomes 



with the small positive correction 5 given by 



Finally, let us consider the orbital system comprising the earth and the moon, for which the 
combined dynamics of the earth-moon system has to be taken into account. 1 This shows [3] 
that measuring the relative earth moon distance by lunar laser ranging, together with the moon's 
orbital period, gives a determination of the combined mass in gravitational units of the earth-moon 
system, which is GM com bined = GM ffl + GM m + GM^ m , with M^ m now the mass of earth-bound 
dark matter 2 lying between the radius of the moon's orbit and the radius of the LAGEOS satellite 
orbit. Thus, by subtracting from GM com bined) as determined by the multi-parameter fit to the 



1 The influence of dark matter on the earth-moon system has been investigated earlier with the aim of determining 
whether the gravitational interaction of galactic dark matter with the ordinary matter in the earth and moon 
obeys the equivalence principle, and the equivalence principle in this context has been verified to high accuracy; 
see [ll| for details. The assumption that dark matter has normal gravitational interactions with ordinary matter 
is implicit in our analysis. 

2 If the earth-bound dark matter distribution is not spherically symmetric, then there will be dark-matter contribu- 
tions to higher multipoles of the earth's gravitational potential, as well as to the monopole mass. 
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lunar laser ranging experiment, the values of CriW© determined by LAGEOS and GM m determined 
by the lunar orbiters, one can get a direct determination of GMd m , 

GM dm = GM comhined - GMjj - GM m , (4) 

subject to our assumption that the moon-bound dark matter lying within the lunar orbiter radius 
can be neglected. Alternatively, if one uses an asteroid determination of GM m , the subtraction to 
be performed is 

GM combmed - GM e - ^® ~ GM dm + GM m S = GM dm + ^GAM ffi (5) 

Since M m /M© ~ 0.0123, © gives 

GM combined - GM e - — ^ ~ GM dm + 0.0123GAM e > GM dm . (6) 

Moreover, if one assumes the earth-bound dark matter AM^ relevant for the asteroid orbit to be 
similar in magnitude to the earth-bound dark M dm matter lying between the moon's orbit and the 
LAGEOS orbit, then ([6]) becomes 

GM combincd - GM e - ^® ~GM dm [l + O(.01)] , (7) 

giving a determination of M dm with a potential one percent accuracy, if the quantities on the left 
hand side of (J7J) were known to sufficient accuracy. In fact, with current data, the errors in the 
left-hand side are more significant than the error arising from the unknown term of order 0.01 on 
the right of ©. 

Proceeding now to a numerical evaluation, 3 the best evaluation of GM§ from LAGEOS data is 
[3| GM ffi = 398600.4415±0.0008km 3 s -2 which when converted to a TDB compatible (Barycentric 



Dynamical Time compatible) figure is|l4] GM e = 398600.4356 ± 0.0008km a s~ 2 . For GM combined , 
the lunar ranging fit EP0 in table 1 of [15f |. in which the sun/ (earth +moon) mass ratio was treated 
as a solution parameter, gives M /M com bi n ed = 328900.5596 ± 0.0011, which converts [ijj] (using 
GM Q = 1.32712440018(8) x 10 n kmV 2 ) to GM combined = 403503. 2357±0.0014km 3 s~ 2 . The lunar 
orbiter measurements reported in [9j] give GM m = 4902.84km 3 s -2 , with an uncertainty of around 
±.05km 3 s~ 2 based on a comparison with alternative determinations, but a much more accurate 
value is obtained from the Eros ranging data of [lol j. which gives R®/ m = 81.300570 ± 0.000005, 



3 The cited papers do not give details about how errors are calculated, but the context suggests that they are 
estimated errors based on residuals to multi-parameter model fits. 
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which corresponds, using the LAGEOS value for GM e , to GM m = 4902.8000 ± 0.0003km V 2 . 
Substituting these numbers into (JH) or ([7]) gives 

GM dm -(403503.2357 ± 0.0014 - 398600.4356 ± 0.0008 - 4902.8000 ± 0.0003)kmV 2 
=(0.0001 ± 0.0016)kmV 2 = (0.3 ± 4) x 10~ 9 GM e , 

(8) 

with the dominant contribution to the error coming from the error in M com bi ne d from the lunar 
laser ranging fit. Thus, current data show that the mass of earth-bound dark matter lying between 
the moon's orbit radius ~ 384, 000 km and the LAGEOS orbit radius 4 ~ 12, 300 km must be less 
than 4 x 10 -9 of the earth's mass, at a 1 a confidence level. 5 As the accuracy of lunar laser ranging 
improves, one can expect this limit on M dm to improve. 

I wish to thank Slava Turyshev for inviting me to speak at the workshop From Quantum to 
Cosmos - III, Airlie, VA that he organized, and for a subsequent email giving me the numbers 
and references used in the numerical evaluation of the preceding paragraph. I also wish to thank 
Peter Bender for a helpful conversation at the Airlie workshop, which was the impetus for this 
investigation, and J. M. Frere for email correspondence. This work was supported in part by the 
Department of Energy under grant no. DE-FG02-90ER40542, and I also wish to acknowledge the 
hospitality of the Aspen Center for Physics. 
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Added Note Gary Gibbons [16( has pointed out that if one assumes that there is no dark 
matter bound to the earth, then the comparison of GM^ as determined by LAGEOS, with that 
determined by lunar ranging, gives a bound on possible non- Newtonian modifications to the grav- 

17], [3], [lj| where the use of satellite 



itational force, and he has alerted me to several references 
orbits to restrict non-Newtonian force models has been discussed. To illustrate with the numbers 
employed above in the dark matter discussion, if one assumes G = Gf ar for the G value relevant 
both for lunar ranging and for the asteroid determination of the earth to moon mass ratio, and 
G = G n ear for the G value relevant for the LAGEOS orbit, and takes M dm = 0, then one has 



4 The LAGEOS orbit is usually described in terms of its altitude of 5, 900 km above the earth's surface, which lies 
about 6, 400 km from the earth's center. 

5 If dark matter gravitationally bound to the earth were assumed uniformly distributed between the moon's or- 
bit and the LAGEOS orbit, this bound, if saturated, would correspond to a dark matter density of order 
10 10 (GeV/c 2 )cm -3 , much higher than the galactic halo density or current limits on the density of dark mat- 
ter gravitationally bound to the solar system. Just based on this, however, one cannot make any statements on 
what should have been seen in dark matter detection experiments, since that would require making assumptions 
about the dark matter density profile around the earth, its mass, and its interaction cross section with ordinary 
matter, all of which enter into determining the experimental sensitivity. No such model-dependent assumptions 
enter the purely gravitational analysis given above. 
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G nC arM e = 398600.4356 ± 0.0008kmV 2 , G far (M e + M m ) = 403503.2357 ± 0.0014kmV 2 , and 
i? e/m = M e /M m = 81.300570 ± 0.000005. When combined these give 

(G ncar - G iax )/G = (0.2 ± 4) x 10~ 9 , (9) 

indicating that G can change by at most ~ 4 x 10~ 9 (the same fractional error that appears in 
([8])) between the radius of the LAGEOS orbit and the radius of the moon's orbit. This is a factor 
of five better than the result given given some time ago by Rapp [18]. However, Turyshev 20J, in 
reviewing fits to lunar ranging, which model the earth-moon distance to 4 mm accuracy, notes that 
"analysis of the LLR data tests the gravitational inverse-square law to 3 x 10~ n of the gravitational 
field strength on scales of the Earth-moon distance". 
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